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ASYMPTOTICS OF ACCESSIBILITY SETS ALONG AN ABNORMAL 

TRAJECTORY 

E. Trelat 1 

Abstract. We describe precisely, under generic conditions, the contact of the accessibility set at 
time T with an abnormal direction, first for a single-input affine control system with constraint on 
the control, and then as an application for a sub-Riemannian system of rank 2. As a consequence we 
obtain in sub-Riemannian geometry a new splitting-up of the sphere near an abnormal minimizer 7 into 
two sectors, bordered by the first Pontryagin's cone along 7, called the L°°-sector and the L 2 -sector. 
Moreover we find again necessary and sufficient conditions of optimality of an abnormal trajectory for 
such systems, for any optimization problem. 

Resume. On decrit precisement, sous des conditions generiques, le contact des ensembles accessibles 
en temps T avec une direction anormale, tout d'abord pour un systeme de controle affine mono- 
entree avec contrainte sur le controle, puis comme application pour un systeme sous-Riemannien de 
rang 2. Comme consequence on obtient en geometrie sous-Riemannienne un nouveau decoupage de 
la sphere au voisinage d'une anormale minimisante 7 en deux secteurs appeles secteur L°° et secteur 
L 2 , delimites par le premier cone de Pontryagin le long de 7. De plus on retrouve des conditions 
necessaires et sufnsantes d'optimalite pour une trajectoire anormale de tels systemes, pour un probleme 
d'optimisation quelconque. 
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I. Introduction 



1.1. Abnormal trajectories 

Consider a control system on W 1 : 



x{t)=f(x(t),u(t)), x(0)=x o (1) 



where / : R™ x R' m — > M™ is smooth, xo £ R", and the set of admissible controls U is made of measurable 
bounded functions u : [0, T(u)] — > Q, C R m . 
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Definition 1.1. Let T > 0. The end-point mapping at time T of system (1) is the mapping 

U — > R" 
Et ■ u ^ x u (T) 

where x u is the trajectory associated to u. 

It is a very classical fact that Et is smooth in the L°° topology, if U C L°°([0, T]). 
Definition 1.2. The timex end- point mapping of system (1) is the mapping 

1+xW — > R™ 
1 (T,u) .— Sr(«) 

Definition 1.3. A control u (or the corresponding trajectory x u ) is said to be abnormal on [0,T] if (T,u) is a 
singular point of the mapping F. 

An equivalent definition may be given using the well-known Pontryagin Maximum Principle (see [21]), which 
parametrizes the previous condition : 

Definition 1.4. A control u on [0,T] (or the corresponding trajectory x u ) is said to be abnormal if there 
exists a non trivial absolutely continuous function p{.) : [0, T] — ► R" x R called adjoint vector, such that the 
trajectory x satifies almost everywhere the system : 

x-— ■ - -— — (2) 

dp ' P dx' du 

where H(x,p, u) = < p, f(x, u) > is the Hamiltonian of the system, and moreover : 

H(x,p,u) = a.c. on[0,T] (3) 

Remark 1.5. If a control u is abnormal on [0, T] then it is abnormal on [0, t] for any t € [0, T]. 

Definition 1.6. An abnormal control is said to be of corank 1 if the adjoint vector p{.) is defined uniquely up 
to a scalar multiplier. 

Definition 1.7. Let u be an abnormal control on [0,T], and x u its associated trajectory. The subspace 
im dE t (u) is called the first Pontryagin 's cone at x^(£). 

An abnormal control is of corank 1 on [0,T] if the subspace im dEx{u) has codimension 1 in R™. 

Remark 1.8. An abnormal control on [0, T] is in particular a singular point of the end-point mapping E T . 
Conversely if u is a singularity of the end-point mapping and if moreover the trajectory x u associated to u is 
almost everywhere tangent to its first Pontryagin's cone then u is abnormal. 



1.2. Accessibility sets 

Definition 1.9. Consider the control system (1), and let T > 0. The accessibility set at time T, denoted by 
Acc(T), is the set of points that can be reached from xq in time T by solutions of system (1), i.e. this is the 
image of the end-point mapping Et- 

Let 7 be a reference trajectory on [0, T], solution of (1), associated to a control u. Our aim is to describe 
Acc(T) near 7(T). If u is not a singular point of the end-point mapping, then obviously Acc(T) is open near 
7(T). The situation when u is a critical point of Et has to be analyzed. 
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1.3. Single-input afRne systems 

Consider a smooth single-input affine system in R n with constraint on the control : 

x(t) = X(x(t)) + u(t)Y(x(t)), \u(t)\ < n 
x(0) = x 

where X, Y are smooth vector fields in M™, and u is a scalar measurable function. Suppose that the trajectory 
7 associated to the control u — is abnormal and of corank 1 on [0, T}. Let p(.) be an adjoint vector associated 
to 7. Then the condition H = is equivalent to : 

<p(t),X(j(t)) >= a.e. 

i.e. X e im dE t (u) along 7. 

The aim of this article is to describe precisely the accessibility set Acc(T) at time T near 7(T). 

The basic object we have to study is the so-called intrinsic second-order derivative of the end-point mapping : 

Definition 1.10. The intrinsic second-order derivative along 7 is the real quadratic form : 

K(v)=p(T).d 2 E T (u).(v,v) 

where v <G ker dEx(u). 

We shall make a spectral analysis of the end-point mapping along such an abnormal trajectory. Using 
the formalism and normal forms of [9] we shall represent the intrinsic second-order derivative by an explicit 
differential operator along the abnormal trajectory 7 ; then a spectral analysis of this operator shall lead to a 
precise description of the boundary of Acc(T) near j(T) (Theorem 2.10). In Section 2.2 we apply this result to 
sub-Riemannian systems of rank 2 (Theorem 2.21), and obtain a new splitting-up of the sub-Riemannian sphere 
near an abnormal minimizer 7 into two sectors, in which the behaviours of minimizing trajectories near 7 are 
topologically different (Section 2.2.2). On the other part this theory on accessibility sets leads to find again 
some well-known results on optimality of abnormal trajectories, that we recall and improve slightly in Section 
3. 

2. ASYMPTOTICS OF THE ACCESSIBILITY SETS 

In this Section we describe precisely the boundary of accessibility sets for a single-input affine system with 
constraint on the input near a reference abnormal trajectory. Then we apply our results to the sub-Riemannian 
case of rank 2 in order to get the contact of the sphere with the abnormal direction. As a consequence we obtain 
a splitting-up of the sphere into two sectors near the abnormal minimizer. 

2.1. Single-input afRne control systems 

Consider a smooth single-input affine control system in MP, n > 3 : 

x(t) = X(x(t)) + u(t)Y(x(t)), x(0) = (5) 

with the constraint on the control 

\u(t)\<r, (6) 

Let Acc n (T) denote the accessibility set at time T for this affine system with constraint n on the control. Let 
7 be a reference trajectory defined on [0, T}. In the sequel we make the following assumptions along 7 : 

(Hq) 7 is injective, associated to u = on [0, T}. 



(4) 
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(i2i) Vi G [0,T] if(i) = Vcct {ad k X.Y{-i{t)) / k e N} (first Pontryagin's cone along 7) has codimension 1, 
and is spanned by the first n — 1 vectors, i.e. : 

K(t) = Vect {ad k X.Y ( 7 (i)) / fc = . . . n - 2} 

(#2) Vi G [0,T] ad 2 y.X( 7 (i)) £ K(t). 

(H 3 ) Vi e [0, T] X( 7 (i)) ^ Vcct {ad fe X.F( 7 (i)) / k = . . . n - 3}. 

(tf 4 ) Vie[o,T] x( 7 (i)) e Jf(t). 

In these conditions 7 is abnormal and of corank 1. Actually assumptions (Hi — H3) are generic, see [10]. 
Moreover we get normal forms in which the intrinsic second-order derivative may be represented by an explicit 
differential operator (see [9] ) , in the following way. 

Lemma 2.1. [9] Under the previous assumptions, the system (X, Y) is in a C° -neighborhood of 7 feedback- 
equivalent to : 

n n— 2 „ n n n— 1 .-, n „ 

/o = ^ + L + 2. "'^•'•"•'•'•'•'7777 + ^ x ^ x ^ + ^ 

i= 1 i)i=2 i=l i=l (7) 

/: 3 



&r„_i 



w/iere o n -i,n-i(*) > ok [0, T], and the 1-jet (resp. 2-jet) of Zi, i — 1 . . . n — 1 (resp. Z n ) along 7 is equal to 0. 
Moreover the feedback (<p, a, /3) satisfies : 

(^a,p).(X,Y) = (f J 1 ) 

where fa = ip*(X + aY), f\ = <p*([3Y), and ip is a germ of diffeomorphism along 7 such that : 
(i) ¥ >(a;i,0,...,0) = (xi,0,...,0) 
(") ^ = ( >--->°>*,0) 

and a, /3 are real functions defined in a neighborhood 0/7 suc/i i/iai /3 does not vanish along 7 and = 0. 
The k-jet is defined as follows : 

" d 

Definition 2.2. Let V = ^ Vfe— — a vector field. Since 7 is given by x\(t) — t and Xi(t) = 0, i = 2 . . . n, the 

fc=i ^ 

+00 

component T4 can be written in a neighborhood of 7 as ^ ipVfe, where joV — V/ 7 , and : 

P =i 

n n 

JiVfe = ^a k i (x 1 )x i , j 2 V k = ^2 b ij( x -t) x i x i' ■ ■ ■ 
Set jiV = > jiVfe— — . Then 7J»V is called i/ie fc-jei ofV along 7. 

fe=l 4=0 

Set xi = i + £. In these conditions, the controllable part of the system is (£,£2, . . . ,x n -\), the reference 
abnormal trajectory is 7(t) = (i, 0, . . . ,0), and the intrinsic second-order derivative Eq(v) along 7 is identified 
to : 

„T n — 1 



7 ° ij=2 
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£1 = £2, • • • ,£n-2 = £n-i,£n-i = v, and &(0) = &(T) = 0, i = 1 . . . , n - 1. 
Integrating by parts, it can be written into two different ways, namely either as a quadratic form in £1 or as a 
quadratic form in £ 2 : 

(1) It is equal to Qi/g^Ci), where : 



„ ± n-2 

Qi(£i) = / «i(£i)dt and gi (&) = V kjgk^ 



with bi-ij-i = a ' J + a '' , and where G\ is the following space corresponding to the kernel ker dEr(O) of 
the first derivative of the end-point mapping : 

Gi = {6 / £i 2( "~ 2)) G L 2 ([0, T]), £«(0) = £«(T) = 0, i = . . . n 2} 
Lemma 2.3. TTie quadratic form Qi is represented on G\ by the operator D\ so that : 

where ( , ) is the usual scalar product in L 2 ([0,T]), and : 

ft - = <8 » 

i=l y i,j=l 

(2) It is equal to Q2/g 2 (.&)> wnere : 

= / 92(6)* and <7 2 (£ 2 ) = V 6 J+ i J+ i4 2) a; 2 j) 
and where G 2 is the space corresponding to the kernel of the first derivative : 

G 2 = {6 /4f ( "~ 3)) eL 2 ([0,T]), £«(0) = £«(T) = 0, i = 0...n-3, and C & dt = 0} 

Jo 

Lemma 2.4. T/ie quadratic form Q 2 is represented on G 2 by the operator D 2 so that : 

Q 2 (6) = (6,^26)^ 

w/iere 

Note that Qi(£) = Q 2 (0 and £>i = ~^ D ^- 

Our aim is to make a spectral analysis of these operators Di,D 2 - Unfortunately the spectrum of Di on d is 
empty. Hence we shall enlarge the Sobolev space Gi so that the spectrum is not trivial and that Representation 
Lemmas 2.3 and 2.4 are still valid. We set : 

• F ! = / ^"~ 2) € L 2 ([0,T]), fW(0) = ^W(T) = o, i = 0. . .n - 3} for the operator D x . Endowed 
with the norm ||£i||fi = 2 '||l 2 j F\ is a Sobolev space. 
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• F 2 = {6 / &~ 3) € L 2 ([0,T]), ^(0) = ^(T) = 0, i = 0. . .n - 4} for D 2 if n > 4 (if n = 3, no 
condition is imposed). Endowed with the norm ||^2||f 2 = 3 '*IIl 2 j -^2 is a Sobolev space. 

Definition 2.5. Let i = 1 or 2. We call T a conjugate time of along 7 is there exists & G Fi\{0} such that 

£(2(n-l-i)) eL 2 ([0)T])and £,.£.=(). 

Lemma 2.6. ('see /or instance [13]) For any f G L 2 ([0,T]), i/T is not a conjugate time, there exists & G F, 
unique such that £j 2 (™ _1_ *)) g L 2 ([0, T]) and = /. Let L denote the operator f £j considered as an 

operator from L 2 ([0,T]) into L 2 ([0,T]) ; it is selfadjoint and compact. 

The following Lemma is an improvement of [9] , where only a non strict inequality is proved : 

Lemma 2.7. Let t c (resp. t cc ) denote the first conjugate time of Qi on F\ (resp. Q 2 on F 2 ). We have : 
< t cc < t c . 

Proof. This proof can be found in [9], where only a non strict inequality is proved. It is included only for 
convenience of the reader. 

Let y G Fx such that y( 2 ("- 2 )) g L 2 ([0,T]), and z = y ; then Q\(y) = Q 2 (z). Let A (resp. 11) denote the 
smallest eigenvalue of Di (resp. D 2 ). From spectral properties of compact selfadjoint operators, we have : 

, . f Qi(y) , . f Q*{z) 

A = mf and fi = mi — - 

(y,y) z ^ ( z , z ) 

By definition : 

Qi(y) Q*(z) 



fi{Slz{s)ds)Ht 

where z G F 2 and / Q T z(t)dt = 0. Hence : 

A>infj rp ® 2{Z) /zGF 2 \{0} 

Moreover we get from the Cauchy-Schwarz inequality : 



f ( [ z(s)ds] dt< [ t f z(s) 2 ds 
Jo \Jo J Jo Jo 



and integrating by parts : 

f T ,t T 2 f T 1 ,T T 2 f T 

/ t / z(s) 2 dsdt < — / z(t) 2 dt - - I t 2 z{t) 2 <-z- z(t) 2 dt 

Jo Jo 2 J 2 Jq 2 J 



Therefore : 

From [22] we know that eigenvalues of D\ (or D 2 ) are continuous and decreasing functions of T. Hence if T = t c 
then A = 0, thus \x < and consequently t cc < t c . □ 

Remark 2.8. If n = 3, we have t cc = +00 provided assumptions (H — H 3 ) are fulfilled on K + . 

Remark 2.9. The notion of conjugate time does not depend on the constraint on the control. It comes from 
the fact that the abnormal reference control belongs to the interior of the domain of constraints. 
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The main result is the following. 

Theorem 2.10. Consider the affine system (5) with the constraint (6), and suppose that assumptions (Hq — H±) 
are fulfilled along the reference abnormal trajectory 7 on [0, T] . Let t cc and t c denote the first conjugate times 
associated to 7. Then : 

(1) There exist coordinates (xi, . . . ,x n ) locally along 7 such that in these coordinates : j(t) — (t,0,... , 0), 
and the first Pontryagin's cone along 7 is : K{t) = Vect {g§^-, ■ • • , dx d ~ }\-y- 

(2) IfT is small enough then for any point (x\, . . . , x n ) of Acc v (T)\{j(T)} close to j(T) we have : x n > 
(see Fig. 1). 





Xi 



abnormal direction 



Figure 1. Shape of Acc v (T), T small 



(3) IfT < t c , then in the plane (xi,x n ), near the point (T, 0), the boundary of Acc v (T) does not depend on 
rj, is a curve of class C 2 tangent to the abnormal direction, and its first term is : 

x n =A T (x 1 -Tf + o((x 1 -T) 2 ) 

The function T 1— > At is continuous and strictly decreasing on [0,t c [. It is positive on [0,i cc [ and 
negative on ]t cc , t c [. 

Moreover, if n depends on x\~T then the result is still valid providing : x\ — T = o{n) as xi — > T. 

(4) IfT>t c then Acd 1 ^) is open near 7(T). 

The evolution in function of T of the intersection of Acc v (T) with the plane (x\,x n ) is represented on Fig. 2. 
The contact with the abnormal direction is of order 2 ; the coefficient At describes the concavity of the curve. 
Beyond t c the accessibility set is open. 

Remark 2.11. The coefficient At can be computed in the following way (see [9], and the Proof just below). 
Actually there exists a function J of class C 2 (n — 2) on [0,T] such that D\J = and satisfying the limit 
conditions : 

Vfce {0,...,n-3} J (fe) (0) = 0, J (fe) (T) = ^ 

Then : 

A T = Q 1 (J) (10) 

Remark 2.12. Let us give a geometric interpretation of the role of the quadratic forms Qi, Q2- 

(1) Role of Q\. On the one part note that G\ is dense in F\ and Qi is continuous on Fx for the norm 
; \\f 1 ; hence the sign of Q\ on G\ is the same as on F\. On the other part, from the definition of the 
first conjugate time t c and the extremal properties of selfadjoint compact operators, we get that Q\ j Fl 
is positive definite if T < t Cl and indefinite if T > t c . Hence the same goes for Qi/q 1 - 
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Figure 2 



Now the interpretation of Qi/Q t is the following : it is an equivalent of the coordinate x n {T) as all 
others coordinates are fixed : X\{T) = T,x 2 {T) = ■ ■ ■ = x n -\(T) = 0. 

As a consequence if T < t c then points (T, 0, . . . , 0, e), where e > are attainable at time T. If 
T > t c , then the same goes with e > or e < 0. 
(2) Role of Q 2 . Contrarily to the previous case, the subspace G 2 is a strict subspace of F 2 . However, set : 

H 2 = {6 / 6 (2( "" 3)) G L 2 ([0,T]),6 W (0) - 6 W (T) = 0,i = 0. ..n- 3} 

Then #2 is dense in F2, and Q2 is continuous on F 2 for the norm || \\f 2 - 

Hence using the same reasoning we obtain the following : if T < t cc then Q 2 / H2 is positive definite ; 
if T > t cc then Q2 /# 2 is indefinite. 

Now note that Q 2 / H2 is an equivalent of the coordinate x n (T) as coordinates x 2 (T), . . . , x n _\{T) are 
fixed to 0, but the coordinate X\{T) is let free. 

As a consequence if T < t cc then points (T ± S, 0, . . . , 0, e), where e, S > are attainable at time T. 
If T > t cc , then the same goes with e > or e < 0. 
This gives us the qualitative shape of Acc{T). Actually this interpretation is enough to study the time- 
optimality of the abnormal trajectory 7, see Section 3.1.1. Here Theorem 2.10 gives more : it describes the 
contact of Acc(T) with the abnormal direction. 

Proof of Theorem 2.10. We shall recall briefly the formalism used by [9], which leads actually to a very precise 
decomposition of the intrinsic second-order derivative. It is based on normal forms of Lemma 2.1. 
Let D\ denote the operator (8) and D 2 the operator (9). 

Lemma 2.13. For all t e [0,T] there exists an Hilbertian basis (e^)„ e N o/L 2 ([0,T]) such that 

• Mn e N e/ n is &{n - 2) and e t n £ F\. 
. Vn £ N D\el - \\e\- 

• Ai < \\ < ■ ■ ■ < \i < ■ ■ ■ 

The eigenvalues \ l n are continuous and decreasing functions oft, and — ► +00. Moreover, if t c denotes 
the first conjugate time of D\ , then : 

• if0<t<t c then A* > 0, 

• ift>t c then X{ < 0. 
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Lemma 2.14. Let < t < t c . Then there exist J*, J* in (^(n — 2)([0, t\), i = 1 . . . n — 2, uniquely defined by 
the following equations : 

(i) Df Jf = D\J\ = 0, i = 1 . . . n - 2 

jW( ) = jW(f) =0 
W j(*)( t ) = ^)(o)=^_ 1) fc = 0...n-3 
Every £ in C^(n - 2) ([(),*]) can 6e expanded in a uniformly convergent series : 

n—2 n — 2 oo 

z— 1 i—1 n—1 

where y n E R, a, - £ (i_1) (0), A = 

We shall use this lemma in the following way. In first approximation the system written in the normal form 
of Lemma 2.1 is given by : 

£ = x 2 , x 2 — x 3 , . . . , x n -i = u 
except the last coordinate, which is given at order 2 by : 



ra-l 

H a 

i,3=2 



ij X j 



Now the meaning of the Ji, Jj's is clear : the functions Ji (resp. Jj) represent the initial (resp. final) conditions 
of the n — 2 first coordinates. More precisely we have, see [9] : 



n-2 



Lemma 2.15. Let T < t c , and £ = ^ VieJ + ^ onJj . We have : 

i=l i=l 

• Efi^-V) = 5\T + a l + odl^-^lk-), i = 1 • • -n - 2. 

n — 2 oo 
i,j=l i—1 

and Afj = Qj(jf , Jj), where Qj denotes the symmetric bilinear form associated to the quadratic form 
Q\. In fact : Qf(£) = d 2 E^{0).(u, u), where u ~ ^ n ~ 1 \ i.e. 

q?(o= r^bijm^m^m (n) 

where fe«-2,n-2 is continuous and strictly positive on [0,T]. 

Hence with this formalism we study the image of the end-point mapping by constructing directly the tra- 
jectories (and not by choosing controls). Actually once a function £ has been determined, the corresponding 
control is u = £, (n ^ + o(\\&^ 1] 

First of all using expression (11) it is easy to prove point 2 of Theorem 2.10. Indeed to study the accessibility 
set at time T from we have to consider functions £ such that : 

£ (i) (0) = 0, i = l...n-2 
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Hence using repeatedly Poincare's inequality 1 we check that if T is small enough then there exists j3 > such 
that : 

QT{O>0 I (&- 2) (t)) 2 dt 

Jo 

We shall deduce the result by proving that actually E^(u) > for the corresponding control u. The result 
would be immediate if E% were CT in H~\ i.e. in L topology for coordinate x n -\. But this is wrong (for 
instance the function / Jo(f 2 + f 3 )dt is not C 2 on L 2 ). Anyway E% is a little better than C 1 . Precisely we 
have (see [4], Theorem 7.1) : 

Lemma 2.16. E%(u) = QT(Q + R(Q where \R(£)\ < C\ \^ n ~ 2 ) | | L » 1 |£(™- 2 ) | \\ 2 . 

And hence if T is small enough then for any point (x\, . . . , x n ) belonging to Acc v (T)\{(T 7 0, . . . , 0)} and to 
the C°-neighborhood of 7 of Lemma 2.1, we have : x n > 0. 

So Lemma 2.15 describes the n — 2 first coordinates and the last coordinate. Anyway [9] do not control the 
coordinate i?J_ 1 (u). Indeed their method consists in replacing the control it e L 2 by the control x n -i G L 2 , 
which is called Goh transformation and consists in fact in immersing L 2 into the space if -1 , see [5]. Such 
trajectories are called generalized trajectories, and this corresponds to add to the set of controls L ([0, T]) Dirac 
measures at and T. 

Here, in order to take into consideration the constraint |u| < 77 and to control the coordinate E n _i(u), the 
method used in [9] has to be adapted. The aim is to describe the boundary of the intersection of Acc v (T) with 
the plane (x\,x n ), close to the point (T, 0) (corresponding to 7(T)). Let x be a real close to T. We have to 
solve equations : 

Ef(u) = x, Ej(u) = ■■■ = El^u) = (12) 

in the domain |u| < 77, and then among such solutions we have to minimize the last coordinate E^(u) (in order 
to describe the boundary of the accessibility set). We proceed in the following way. Set x\{t) = t + £(t) ; we 
shall first determine £ such that the control u = £( n ~ 1 '> satisfies (12). From Lemma 2.14, £ can be expanded in : 

71 — 2 00 

£ = ^a;jf + ^>ief (13) 

i=l i=l 

n-2 

From Lemma (2.15), the first n — 2 coordinates are, if u ~ £( n_1 ) = \^ a^jf^ ^ + v : 

i=l 

Ef(u) = T + ai+o(ai,...,a n _ 2 ,|H|L~) 
Sj(tt) = a 2 + o(ai,...,a n _2,||u||L«>) 

E n-2( u ) = an-2+o(ai,...,a n _ 2 ,||u||L~) 

Let us solve equations (12). We get from the Implicit Function Theorem : 

ai - 2;-T + o(a;-T,||v|| L ~) 
a 2 = o(s-T,||u|| L oo) 

OLn-2 = 0(X - T,\\v\\ L co) 



1 /(0) = =^ ||/||l~ ( [o,t]) < V5 5 ||/ , || L "([o,ri) 
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as x — > T and v ^> 0. Further we will modify u in order to get also : 

Zn-i(0) = x n -i(T) = and |u| < 77 

Let us first study the last coordinate. We get from Lemma (2.15) : 

00 

El{u) =A T (x- T) 2 + ]T Xfvf + o((x-T+\ \v\\^ ) 2 ) 

i=l 

as x — > T and i> > 0. In this expression we can see that the minimum is reached at v — up to a o((x — T) 2 ), 
for if T < t c then all Xf 's are positive (Lemma 2.13). 

Remark 2.17. If we do not neglect any term then the minimum is not necessarily reached at v = because of 
terms of order 3 in (x — T) 2 Vi. 

n-l 

Therefore, at order 2, the minimum is reached at u <~ ^ ctijj and equals : 

i=l 

El{u) = A T {x - T) 2 + o((x - T) 2 ) asi^T 

We shall now prove that it is possible to modify u, taking into account the constraint \u\ < 77, in order to 
obtain x n -i(0) = x n -\{T) — 0, without changing the previous results on the other coordinates. 
First of all, if x is close enough to T, then the aj's are small and thus the constraint \u\ < r\ is satisfied. 
Let us modify u (and hence x n -i) in the following way. Set u(t) = 77 if x n -i(0) > or —rj if x„_i(0) < on 
[0,ti] (same construction on [t2,T]), where t\ and t 2 are the coinciding times, i.e. the times at which x n -\ 
coincides with its initial graph, see Fig. 3. 



slope rj 




-slope — 77 



Figure 3 



Precisely : 



n-2 



if < t < ti : \u(t)\ = r), |a; n _i(i)| = r)t, and the coinciding time t\ is such that 77*1 ~ | ^^aijf^ n 2 \t\ 



n-2 



if ti < * < *2 : u{t)~^-V{t),x n - 1 {t) = Y,UiJl {n 2 \t)- 



n-2 



if t 2 < t < T : \u(t)\ = 7/, |x„_i(t)| = t/(T - t), where t 2 is such that t?(T - t 2 ) ~ | ^ «iJf ( " 2) (* 2 ) | 
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We shall now check that if 

x — T = o{n) asx^T 

then all previous results are still valid. 

Clearly : nt\ — 0(x — T) and n(t2 — T) = 0(x — T) as x — > T. Moreover, at order 1, the system is : 

±i = l + x 2 , x 2 = x 3 , . . . , x n -i = u 

Therefore : 

.71— 1 .n— 1 



-x, and a;i(i 2 )-i 2 -a;i(T)+T~r ? - 



>-l)! 1V ^ 7 ' 1V 7 (n-1)! 

Now if n > 3 and x — T = 0(77) then these terms are negligibly small in comparison to x — T. 
As concerns the last coordinate, we obtain : 

x n (h) = 0(^4) = o((x - T) 2 ) and x n (T) — x n (t 2 ) = 0(n 2 tl) = o((x - T) 2 ) 

In these conditions, all our previous construction is still valid. Hence in the plane (xi,x n ) the boundary of 
Acc T, (T) is a curve of class C 2 , independant of the constraint, such that x n ~ At{x\ — T) 2 , which proves 
Theorem 2.10. Moreover it results from [9] and [22] that the function T At is continuous and decreasing on 
[0, t c [, positive on [0, t cc [ and negative on ]t cc , t c [. □ 

2.2. Application to the sub-Riemannian case 

2.2.1. Asymptotics of the sub-Riemannian sphere along an abnormal direction 

Consider a smooth sub-Riemannian structure (M, A, g) where M is a Riemannian n-dimensional manifold, 
n > 3, A is a rank 2 distribution on M, and g is a metric on A. Let xq G M ; our point of view is local and we 
can assume that M = R" and x = 0. Suppose there exists a smooth injective abnormal trajectory 7 passing 
through 0. Up to changing coordinates and reparametrizing we can assume that : 

. 7 (t) = (i,0,...,0), 

• A = Span {X, Y} where X, Y are g-orthonormal, 

• 7 is the integral curve of X passing through 0. 

Under these assumptions, the sub-Ricmannian problem is equivalent to the time-optimal problem for the system : 

x = vX(x) + uY(x), x(0) = (14) 
where the controls v, u satisfy the constraint : 

v 2 +u 2 <l (15) 

The reference abnormal trajectory 7 corresponds to the control : v = l,u = 0. 
Let us now define a notion of constrained accessibility set : 

Definition 2.18. Let < a < 1. We denote by Accg R (T) the accessibility set at time T for the sub-Riemannian 
system (14) with the additional constraint on the control : 



v 2 + u 2 < 1 , 1 — a < v < 1 , lul < a 



(see Fig. 4) 



Note that controls steering to points of Accg R (T) are in a a-neighborhood in L°° metric of the abnormal 
reference control v — 1, u = 0. 
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Figure 4 

Definition 2.19. We call affine system associated to the sub-Riemannian system (14) the following system : 

x = X{x) +wY(w) (16) 

where the control w satisfies a constraint of the form : |tu| < t]. 

Let Acc v A (T) denote the accessibility set at time T for this affine system with the constraint : \w\ < rj. The 
reference trajectory 7 corresponds to w = 0, and is also abnormal for this affine system. 

The following lemma gives a precise comparison of constrained accessibility sets of systems (14) and (16) : 

Lemma 2.20. (l)Vae]0,l[ Acc a SR {T) <Z (J Accf^ ' (s) 

(l-a)T<s<T 

To V \ 



(2) 

T„<s< 



(J Acc^is) C Acc% R (T) where a = max 1 



Proof. Let us prove the first inclusion. If x\ G Accg R (T) then there exists a control (v, u) such that v 2 + u 2 < 1, 
1 — a < v < 1, |u| < a, and such that the corresponding trajectory satisfies : 

x = vfo + ufi, x(Q) = 0, x(T) = xi 

As a < 1, v does not vanish, the following reparametrizing holds : ^| = v. Set y(s) — x(t), w(s) = -^y, and 
S = J T v. Then : x\ = x(T) = y(S), and S is such that : 

(1 - a)T < S < T 

Moreover : 

d y t , u r t j f 
-j" = Jo + - Ji = Jo + w/i 

as v 

where \w\ < -r 2 — • Therefore : 

I I — 1 — a 

Xi e (J ^cc^(s) 

(1-q)T<;s<T 

which proves the first part of the lemma. 

Let us now check the second inclusion. Let S G [To, . T ] and xi G Acc^S*). There exists a control w 

such that \w\ < rj and the corresponding trajectory satisfies : 

dy 



(lg =fo+wf u y(0) = 0, y(S) = Xl 
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Let £ > such that T — Let us make the reparametrizing : ^| = 1 — s, and set : x(t) = y(s), v(t) = 1 — e, 
u(t) = (1 - s)w(s), where t G [0,T]. Then : 

x = vf + u/i, &(0) = 0, ar(T) = y(S) = Xl 

Let us now check the constraint on the control (v,u). By definition : T < S < -y= = and 5 = (1 — e)T. 



Hence : 

I 



1 — e < — , and 1 — e > 5r 



and thus 



1 - I 1 - 1 < v < 1 and Id < 



2 



Moreover : 

v 2 + u 2 < (1 - e) 2 (l + r? 2 ) < 1 

Therefore we can conclude that : 

Xl &Acc% R {T) where a = max ^1 - -j==^ 

□ 

Using the previous lemma and Theorem 2.10 we can prove the following : 

Theorem 2.21. Suppose assumptions (Hq — H 3 ) are fulfilled along the reference abnormal trajectory 7 for the 
system (X, Y). Let t cc and t c denote the first conjugate times of '7 for the associated affine system. Let a e]0, 1[. 
Then : 

(1) There exist coordinates (xi, . . . , x n ) locally along 7 such that in these coordinates : "f(t) = [t, 0, . . . , 0), 
and the first Pontryagin's cone along 7 is : K{t) = Vect {gf^-, ■ • • , dx d ~ }\-y- 

(2) If T is small enough then for any point (xi, . . . ,x n ) of Accg R (T) close to 7(T) we have x n > (see 
Fig. 1). 

(3) If T < t cc , then in the plane (xi,x n ), close to the point (T, 0), the boundary of AcCg R (T) does not 
depend on a, is a curve of class C 2 outside (T, 0), tangent to the abnormal direction, whose first term 
is : 

• if x\ <T then x n = 0. 

• ifx!>T then x n = A T ( Xl - T) 2 + o((xi - T) 2 ). 
The function T At is the same as in Theorem 2.10. 

(4) IfT> t cc then AcCg R (T) is open near 7(T). 

Figure 5 represents the evolution of Accg R (T) in function of T in the plane (xi,x n ). It is open in a neigh- 
borhood of j(T) if T > t cc , contrarily to the affine case where it becomes open only beyond t c . 

Remark 2.22. To compare the system (14) with its associated affine system (16) we need the following 
reparametrizing (see Proof of Lemma 2.20) : 

ds 

~r = v 
dt 

which only holds if v does not vanish. This condition is satisfied when the control (v, u) is in a a-neighborhood 
in L°° metric of the abnormal reference control (1,0), for in this case v is close to 1 in L°°. Hence using this 
method it is only possible to describe a constrained accessibility set, i.e. in a a-neighborhood in L°° metric of 
the reference abnormal control. 
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Figure 5 

Proof of Theorem 2.21. The aim is to compare precisely systems (14) and (16) using Lemma 2.20 and to apply 
Theorem 2.10. In order to do this we first have to normalize the affine system (16) using Lemma 2.1. We denote 
by Acc A (T) the accessibility set at time T with constraint (3 of the affine system (16). 

The system (7), x = fa + ufi, is called normalized affine system and will be refered as (AN). Let Acc AN (T) 
denote the accessibility set at time T for this system with the constraint : |u| < 77. Due to the particular forms 
of the feedback (ip, a, (i) and of the system, we have : 

Lemma 2.23. ip~ 1 (Acc n AN (T)) = Acc A (T) where (3 = 0(r]) asn^Q. 

We know from Lemma 2.20 that : 

Acc a SR (T)c |J Accf^(s) 

s<T 

Hence in the normalized coordinates we get : 

<p(Acc a SR (T)) c |J Acc p AN {s) (17) 

s<T 

We shall use Theorem 2.10, which describes the boundary of Acc AN (s) in the plane (xi,x n ), to study the 
boundary of Acc AN (s). Using the fact that t The function 1 1— > At is continuous and decreasing on [0, t cc [, 

s<T 

we can assert that the boundary of Acc AN (s) is given, in the plane (x\,x n ), close to the point (T, 0), by 

s<T 

the following curve, see Fig. 6 : 

• if x\ < T then x n = 0. 

• if X! > T then x n = A T (xi - T) 2 + o((a;i - T) 2 ) 

Let x n — f(x\) denote the curve parametrizing the boundary of <p(Acdg R (T)) in the plane [x\,x n ). From 
inclusion (17) we get : 

• if x x < T then f(x 1 ) > 0, 

• if xi > T then f( Xl ) > A T (xi - T) 2 + o((xi - T) 2 ). 

Let us now prove the converse inequality. To this aim we shall use varying constraints depending on x± — T . 
We proceed in the following way. From Lemma 2.23 : 

y- l (Ace AN (^L=)) c Acc A ( - T ) 
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Figure 6 



where (3 = 0(rj). Without loss of generality we can assume (5 < rj, and thus : 

Ac^ ( -?=L=) C Acc\ ( -?=L=) 



Now from Lemma 2.20, we get that for any rj > small enough : 

T 



^( 7 _)CA CC ^(T) 



And thus in the normalized coordinates 



Acc\J^==) C V {Acc% R {T)) (18) 
Let /an(xi) be the function parametrizing the boundary of Acc\ N (^ -j= = ^ in the plane (x\,x n ). We know 



from Theorem 2.10 that : 



Jan{xi) = A T (xi - 



i + v 2j W^/T^fJ ) 



provided x\ 



T 



V 



0(77). This latter condition is fulfilled if r] = (xi — T)* , and in this case we have 



y/ 1 + 7 f 

moreover : rj 2 = o(x\ — T). On the other part, from the continuity of t i— > A t : 

A T = A T + o(l) asxi->T 

We obtain : Jan(xi) = At(xi — T) 2 + o((x\ — T) 2 ). Finally, from inclusion (18) we conclude : 

f(xi) < A T ( Xl - T) 2 + o((xi - T) 2 ) 

which ends the proof. 



□ 
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2.2.2. Splitting-up of the sphere near an abnormal direction 

Let T > small enough so that properties 2 and 3 of Theorem 2.21 are satisfied. In particular the reference 
abnormal trajectory 7 is minimizing, see Section 3.1. Then A = j(T) belongs to the sub-Riemannian sphere 
5(0, T) with radius T. If controls steering to points of the boundary of AcCg R (T) in x n > (that are L°° - 
optimal) are actually globally optimal, then this boundary is included in the sphere 5(0, T). In this case the 
sphere splits into two sectors near j(T), bordered by the first Pontryagin's cone x n = : 

• sector x n > corresponding to the previous description, 

• sector x n < 0. 

According to the previous results, final points at time T associated to controls which are L°°-close to the 
reference abnormal control are in the first sector : x n > 0. Obviously due to controllability of the system the 
sector x n < is accessible. In fact a basic calculus shows : 

Lemma 2.24. For any neighborhood V of the point A in M™ we have : 

s(o,T)nvn(x n < 0) ^0 

These points in [x n < 0) are reached by controls which are close to the reference control in L 2 metric but 
not in L°° metric. More precisely : 

Lemma 2.25. Let M n = Et{u 71 ) € 5(0, T) whose last coordinate x n is strictly negative. Let u denote the 
abnormal reference control. We suppose that M n converges to A = E{u) in R". Then u n converges to u in 
L 2 ([0,T]) but not in L°°([Q,T}). 

Hence near the abnormal direction the sphere is splits into two sectors : 

• the L°° -sector : (x n > 0) fl 5(0, T) (described by Theorem 2.21), made of end-points of minimizing 
trajectories associated to controls L°°-close to the abnormal reference control, 

• the L 2 -sector : (x n < 0) n 5(0, T), made of points reached by minimizing controls L 2 -close, but not 
L°°-close to the abnormal reference control. 

The contact of the first sector is known, but not the second one a priori. Anyway according to the Tangency The- 
orem (see [23]), under some nice stratification assumptions, this L 2 -sector ramifies tangently to the Pontryagin 
cone x n = 0, see Fig. 7. 




Figure 7 
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Remark 2.26. In particular, minimizing trajectories 700 joining to points of the sphere in the L°°-sector are 
C 1 -close to the reference abnormal trajectory 7. Minimizing trajectories 72 joining to points of the sphere in 
the L 2 -sector are C°-close, but not C 1 -close to 7, see Fig. 7. 

Typical example : the Martinet case. Consider the two following vector fields in R 3 : 

X = — + ,y = — 

dx 2 dz dy 

and endow the distribution spanned by these vector fields with an analytic metric g of the type : 

g = adx 2 + cdy 2 

where a = (1 + cay) 2 and c = (1 + f3x + 7y) 2 . The abnormal reference control for the sub-Riemannian system 
x = vX(x) + uY(x) with constraint v 2 + u 2 < 1 is v = l,u = 0, and corresponds to the trajectory 7 : 
x(t) = t, y(t) = z(t) = 0. We have, see [11] : 

Lemma 2.27. Assumptions (Hq — H3) are fulfilled along 7 if and only j/a^fl. In this case branches 1 and 2 
(see Fig. 7 with x\ — x, x n = z) have the following contacts with the abnormal direction : 

• branch 1 : x > T, z= ^i{x - T) 2 + o((x - T) 2 ) 

• branch 2 ; x < T, z ~ |(1 + 0(T))(x - T) 3 

Remark 2.28. The coefficient At of the first branch can be computed directly or using formula (10) (see 
Remark 2.11). 

As we are in dimension 3, results of Theorem 2.21 are in fact available on M+, see Remark 2.8. The L 2 -sector 
is z < and the L°°-sector is z > 0. 

Proof of Lemma 2.24- Consider the sub-Riemannian system (14) in the coordinates (x\, . . . ,x n ) of Theorem 
2.21 : 

x 1 = v(l + x 2 + x n f 1 (x 1 ) + Zi) 
X2 = v(x 3 + x n f 2 (x 1 ) + Z 2 ) 

X n -\ = U + X n f n -i(xi) + Z n _i 
n 

X n — V ^ ^ 0{j {x^jx^x j -\- Z n 
i=2 

where the 1-jet (resp. the 2-jet) of Z\, . . . , Z n _i (rcsp. Z n ) along u — is equal to 0. 

Our goal is to construct a control close to the reference abnormal control in L 2 metric, whose associated final 
point at time T is close to (T, 0, . . . , 0) and is such that x n < 0. Let e > 0. Consider the following control v 
(see Fig. 8) : 

• if < t < ^ then v(t) = 1, 

• if < t < Z±£ then v(t) = -1, 

• if ^ < t < T then v(t) = 1. 

Set u = 0. It is clear that (v,u) is abnormal (but not minimizing). Consider the following perturbation (see 
Fig. 9) : 

• if < t < ^ then Sv(t) = Su(t) = 0, 

• if < t < I then 5v(t) = 1 - VT^s 2 , 5u(t) = e, 

• if f < t < Z±z then Sv(t) = 1 - VT^e 2 , 5u{t) = -e, 

• if 2i±i <t<T then Sv{t) = Su(t) = 0. 
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v(t) 



Figure 8 



1 - VT^e 1 



Sv(t) 



T 



6u(t) 



Figure 9 



It is clear that (v + Sv) 2 + (6u) 2 = 1. Moreover : \\5v\\ L i = 0(e 3 ) and ||5u|| L i = 0(e 2 ). 
The end-point mapping Et being C°° in L 1 topology (see for instance [17]), we have : 

Vi e {1, . . . ,n - 1} E l T (v + Sv, Su) = E l T {v, 0) + dE l T (v, 0).(6v, Su) + 0(e 4 ) 



and 



E%(v + Sv, Su) = 0) + dE%(v, 0).(Sv, Su) + ^d 2 E%(v, 0).(Sv, Su)® + 0(e 6 ) 



Moreover we have Ej,(v, 0) = 5}(T — e) for i = 1, . . . , n. On the other part the linearized system along (v 

yi=Sv + vy 2 , V2 = vy 3 , y n - 2 = vy n -\, y n -i = Su, y n = 
Let us calculate y n -i(t) '■ 

2 



if < t < ^ then y n -i(t) = 0, 



• if < t < f then y n -i(t) = e(t - 

• if f < t < ^ then y n -i(t) = -s(t - 2^£), 

• if < t < T then y n -i(t) = 0. 

In particular : ||t/ n -i||L~ = 0(e 2 ), and y n -i is equal to outside an interval of length e. We get easily 



|2M-2||l= = • • • = = 0(e 3 ) 
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Hence using (19) : 

Ef(v + Sv,Su) = T- £ + 0(e 3 ) 
Ef(v + Sv,Su) = 0(e 3 ) 

El_ 2 {v + 5v,5u) = 0(e 3 ) 
*%_i{v + 5v,5v) = 0(e 2 ) 
Let us now compute the last coordinate. We have to calculate z n , where : 



-l 

E 



v V a l0 (t)y t yj 



It represents the intrinsic second-order derivative. From (21) we get : 

z n (t) = Hi)«„-i,„-i(i)!/li(i) + y n -i(t)0(e 3 ) + 0(e 6 ) 
As y n -i is equal to outside [^-j 2 , ^J 2 ], we have : 



/ 

Jo 



T 

3 



y n _i(t)0(e d )dt - 0(e b ) 



Moreover the coefBcient a„_i jn _i is continuous and does not vanish on [0, T], hence there exists a > such 
that a„_i in _i(t) > a on [0,T]. We get : 



z n (T) < -a ^ i£_!(t)dt + 0( £ 6 ) < + 0( £ 6 ) 

Therefore from (20) : 

Hence for any neighborhood V of (T, 0, . . . , 0), j4cc(T)nF contains points such that x n < 0, and hence the same 
goes for S(0, T) n V since the abnormal reference trajectory is minimizing. By construction, controls steering 
to these points are L 2 -close but not L°°-closc to the reference abnormal control. □ 

Proof of Lemma 2.25. As M n belongs to the sphere 5(0, T), we have : ||u„|| L 2 = ||u||l 2 = T, thus the sequence 
(u n ) ne N is bounded in L 2 . Hence up to a subsequence we can assume that u n converges weakly to v G L 2 
(denoted by u n <^>v). From the continuity of the end-point mapping Et in the weak topology on L 2 (see [23]), 
we can assert that M n — E T{u n ) converges to Et{v), and thus Et{v) = Et{u). The assumptions on the 
reference abnormal trajectory imply that v = u. Hence u n ^> u, and on the other part : ||u n ||L 2 = II u IIl 2 ; 
therefore u n converges (strongly) towards u in L 2 . 

Moreover results stated by Theorem 2.21 imply that u n docs not converge towards u in L°° (because it would 
imply that x n > 0). □ 



3. Application : optimality of abnormal trajectories 

In this Section we apply our previous theory on accessibility sets to studying optimality of abnormal tra- 
jectories ; this leads us to find again some well-known results. Indeed in the notations of Theorem 2.10, this 
Theorem implies in particular that if T < t cc then 7 is isolated in CP -topology in the space of all trajectories 
which connect given end-points, and thus is optimal for any cost in this topology. This well-known property, 
called rigidity, was intensively studied. The main results concerning this analysis in a generic context were given 
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first in [22] and [9] for single-input affine systems, then in [4,14,18,25], for sub-Riemannian systems, and in [6] 
in general. Moreover these authors developed a Morse theory in order to characterize conjugate points, that is, 
points beyond which the abnormal trajectory is no more optimal. 

Hence results given in this Section are not really new. However they are slightly different from the results cited 
above. Indeed on the one part in [9] were obtained necessary and sufficient conditions for C°-time-optimality 
of abnormal trajectories of single-input affine systems without any constraint on the control. Here we improve 
their statement by adding a constraint on the control and studying the problem of minimizing any cost. On the 
other part, in [6] was made a general theory (i.e. for nonlinear systems) on optimality of abnormal trajectories 
in L°° topology on the controls. Results given here are valid in the C°-topology on the trajectories (but only 
for single- input affine systems). Moreover we study the equivalence between the time-optimality problem and 
the problem of minimizing any cost, the final time being fixed or not. Finally Theorem 3.10, which concerns 
optimality of abnormal trajectories for sub-Riemannian systems of rank 2, makes a link between the works of [9] 
and [4,5]. 

3.1. Optimality of abnormal trajectories for single-input affine systems 

Consider the single-input affine system (5) with constraint (6), and suppose assumptions (Ha—Hi) are fulfilled 
along a reference abnormal trajectory 7. We first study the time-optimal problem, and then the problem of 
minimizing some cost. 

3.1.1. Time optimality 

Definition 3.1. • The trajectory 7 is said C° -time-minimal on [0, T] if there exists a C°-neighborhood 

of 7 such that T is the minimal time to steer 7(0) to j(T) among the solutions of the system (5) with 
the constraint (6) that are entirely contained in this neighborhood. 
• Recall that 7 is associated to the control u = 0. Let S > 0. The trajectory 7 is said L°° -time-minimal 
on [0, T] if there exists a neighborhood of in L°°([0, T + 5]) such that T is the minimal time to steer 
7(0) to 7(T) among trajectories associated to controls of this neighborhood. 
Obviously if 7 is C°-time-minimal then it is L°°-timc-minimal. 

Theorem 3.2. Under assumptions of Theorem 2.10, the trajectory 7 is Cf -time-minimal if and only ifT < t cc . 
Moreover 7 is not L°° -time-minimal if T > t cc . 

The proof is clear by inspecting Fig. 2 and Proof of Theorem 2.10. 

Remark 3.3. If n — 3, we have t cc = +00 provided assumptions (Ho — H3) are fulfilled on R + . Hence in this 
case 7 is C°-time- minimal on R + . 

3.1.2. Optimality for some cost 

Let us now consider the problem of minimizing some cost C(T,u), also denoted by Ct(u), where C is a 
smooth function satisfying the following additional assumption along the reference singular trajectory 7 : 

(H 5 ) VT rank (dE T (0),dC T (0)) = n 
i.e. the singularity of the end-point mapping of the extended system has codimension 1, and in particular the 
cost is independant from the end-point mapping along 7. We consider several optimization problems : 

(1) final time not fixed : the aim is to steer the system from xo to x\ in some time T (not preassigned) and 
minimizing the cost C. 

(2) final time fixed : let T > fixed ; the aim is to steer the system from x$ to x\ in time T and minimizing 
the cost Ct- 

1. Final time not fixed. 

Definition 3.4. • The trajectory 7 is said to be C° -cost-minimal on [0, T] if there exists a C°-neighborhood 

of 7 such that for any trajectory q contained in this neighborhood, with q(0) = 7(0) and q(t) — j(T), 
we have : C(t,v) > C(T, 0), where v is the control associated to q. 
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• Let S > 0. The trajectory 7 is said to be L°° -cost-minimal on [0, T] if there exists a neighborhood of 
in L°°([0, T + 5}) such that, for any trajectory q associated to a control v of this neighborhood, with 
q(0) = 7(0) and q(t) = 7 (T), we have : C(t,v) > C(T,0). 
Obviously the C°-cost-minimality implies the L^-cost-minimality. 

We have the following result (compare with [6]) : 

Theorem 3.5. Under assumptions (Hq — H5), the trajectory 7 is C° -cost-minimal if and only if it is C° -time- 
minimal. Actually, 7 is C° -cost-minimal ifT< t cc , and is not L°° -cost- minimal ifT> t cc . 

Hence if the final time is not fixed then both problems of cost-minimization and time-minimization are 
equivalent. 

Proof. If T < t cc , then the abnormal trajectory 7 is isolated in a C°-neighborhood, hence in particular is C°- 
cost-minimal. If T > t cc , we know that 7 is no more time-minimal since the point (T, 0, . . . , 0) belongs to Acc(t) 
for some t <T, see fig. 10. 




Figure 10 

To prove that 7 is no more L°°-cost- minimal, we have to show that the point (T, 0, . . . , 0) belongs to M r {t) 
for some r < Ct(0), where Ct(0) is the cost of 7 at time T, and M r (t) is the level set r at time t of the value 
function associated to the cost C. 

In the coordinates (x\, . . . , x n ) of Lemma 2.1, the end-point mapping is E t = (E{, . . . , E^), where : 

rank {dE{{Q), . . . , dE f n (0)) = n - 1 and dE^ (0) = 

Instead of working in L 2 ([0, T]), we choose as control space the Sobolev space H 1 ([0, T]) = {/ e L 2 ([0, T]) / /' e 

L 2 ([0,T])}. Endowed with the norm H/Hh 1 — v/II/IIl 2 ~^ II/'IIl 2 ' this is a Hilbert space. The reason to use it 

is the following : it can be immersed in a compact way in C°([0, T]), and in the proof we shall need continuous 
controls. 

The end-point mapping is still differentiable in H 1 since it is differentiable in L 2 . From assumption (H4), 
the end-point mapping is independant from the cost along 7, hence at time T there exist n independant vector 
fields , ef , . . . , e^^i such that : 

dC T {0).e% = 1 and dEf(0).eJ = 1, i=\...n-\ 

We can decompose H 1 in the following way : 

H 1 ([0, T]) = Ref © • • • © Re^ © ReJ © F 
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where ker dE? (0) = ©Me J © F and ker dC T (0) = 0KeJ F. 

Let Qt(u) = d 2 E T (0).(u, u) (as in Lemma 2.15). The kernel of the first derivative is : ker dEr(0) = KeJ ©F, 
and hence the intrinsic second-order derivative at time T along 7 is the quadratic form Qt j R e T eF- O n the 
other part the space Kef ©Reg 0F represents the domain of the latter quadratic form where the limit condition 
on the first coordinate has been relaxed, hence it is the domain of the reduced operator D\. Now by definition 
of the first conjugate times t cc and t c , we have : 

• if T < t cc : Q T/ MeJme^eF is positive definite, 

• if t cc < T < t c : Q T 1 R e r ffiF is positive definite, and Qt / R e T ffiRe T ffif is indefinite, 

• if T > t c : Qt 1 m.^®f is indefinite. 

From now on we assume that T > t cc . Then there exist v,w in ReJ © F such that Qt{z[ + v) < and 
Qt(w) > 0. 

Remark 3.6. Note that we do not have necessarily Qt{^[) < 0. 

Let S > small. As v and w are continuous, it is possible to extend them on [T, T + 5] respectively by v(T) 
and w(T) (here is the role of H 1 ). Consider now the following control u € H 1 ([0, T]) : 

u = a\e\ + a 2 e^ H V a n -\e^_ x + a\v + a n w 

where a\, . . . , a n are real numbers. In the same way we can extend u on [0, T + 5}. We shall prove that we can 
choose ai, . . . , a n small and t close to T such that the previous control u satisfies : 

M < V , t<T + S , E t (u) = 7(F) , C t {u) < C T (0) 

i.e. 7(F) belongs to M r (t) for some r < Ct(0). To this aim we shall first use the Implicit Function Theorem on 
a\, . . . , a n (to take into account the n — 1 first coordinates, i.e. the controllable part of the system), then the 
Mean Value Theorem on the last coordinate, and endly choose t in order to make the cost lower than Ct(0). 
Note that in this method all reasonings are in finite dimension. 
Let us expand, for t close to F and i = 1, . . . , n — 1 : 

El(u) = F<(0) + dEl(0).u + o(||«||Hi( [0 , t ])) 

Moreover : 

• E\ (0) = t, and for i = 2, . . . , n - 1 : F*(0) = 0, 

• o(||u|| H i([o,t])) = o(II«IIh 1 ([o,t+ < 5])) = o(ai, . . . ,a„), 

• For i = 1, . . . ,n - 1 : dEj(0) = dE?{0) + 0{t - T). 

Therefore we get, for i = 1, . . . , n — 1 : 

dE\(Q).u = a t + o(ai, . . . , a n , t - T) 

and thus : 

E t (u) = (t + ai+ o(ai, ...,a n ,t — F), a 2 + o(ai, ...,a n ,t — F), 

a n _i+o(ai,...,a n ,t-T), F*(u)) 

Let us first solve the system of equations : 

E[{u)=T, E t 2 (u) = --- = E t n _ 1 (u) = 
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We get from the Implicit Function Theorem : 

V*e{l,...,n-1} Oi = fi(T-t,a n ) 

and moreover : 

a\ = T — t + o(T — t, a n ) and for i = 2, . . . , n — 1 : dj = o(T — a n ) 
Let us now calculate the last coordinate : 

Ei(u) = ££(0) + d££(0).u+ld 2 ££(0).(u,u) +o(||«||| 1 ) = Q t («) + o((T - t) 2 , a 2 ) 

Moreover : 

<9t( u ) = Qt(ai(ei +v) + a 2 e^ H h a n -\e^_ x + a n w) 

Hence 

££(u) = (T — t) 2 Q t (ef + «) + a 2 n Q t (w) + a n (T - t)Q t (eJ + v,w)+ o((T - i) 2 , a 2 ) 

= g(t,a n ) 

where Q denotes the symmetric bilinear form associated to the quadratic form Q. On the one part : g(T, a n ) = 
a 2 Qt{w) + °( a n) > if o„ is small enough (by definition of w). On the other part : g(t, 0) = (T — i) 2 Q t (ef + 
v) + o((T — t) 2 ). From the definition of v, we have : Qr{eJ + v) < 0. Hence if t is close enough to T, but not 
equal to T, we get : g(t,0) < 0. Now applying the Mean Value Theorem to g we can assert that there exist t 
close to T and a n small such that g(t, a n ) = 0, that is E* n (u) — 0. Actually there exist four such couples (t, a n ) 
with t > or t < 0, and a n > or o„ < 0. 

Remark 3.7. For such a couple (i, a„) we can say more about the asymptotics of t with respect to a n , which 
will be useful in the sequel. To this aim let us solve the following equation, where a n is fixed : 

g(t, an) = (22) 

We have : 

• Qt(eJ + v) = Q T {e\ + v) + 0(T - t), 

• Qt(w) = Qt{w) + 0(T - t), 

• Q t (eJ + v, w) = Q T (eJ + v,w) + 0(T - t). 

Hence equation (22) becomes : 

Q T {el + v)(T- tf + a n Q T {el + v,w){T -t) + a 2 n Q T {w) + o((T - t) 2 , a 2 ) = 

If a n is fixed, close to 0, we easily find two solutions t\ < T and t^> T such that both satisfy : 

T - t = /3a n + o(a„), G M\{0} 

To end the proof we have to study the cost Ct (u) and prove that t and a n can be chosen so that this cost is 
lower than Ct(0). Recall vectors v, w e Me^ © F previously defined so that : 

Q T (ei+v)<0 and Q T (w) > (23) 



We can write : 



v = XeQ +vi, w = [le^ + w\ where v\,wi G F 
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Up to adding eej, e small, we can assume that fj, ^ 0. This does not change anything in inequalities (23), nor 
in our previous reasoning. Let us denote C(t, u) = C t (u). We can expand : 

8(1 

C t (u) = C r (0) + a(t — T) + — (T, 0).u + o(T - t, \\u\\) 

where a — ^jr(T, 0). Moreover by definition of u and ef : 

dC 

— (T,0).u = dC T (0).u = ai\ + a n n 
ou 

On the one part : a\ = T — t + o(T — t, a n ). On the other part from Remark 3.7 there exists [3 ^ such that 
a„ = /3(t — T) + o(t — T). Hence we obtain : 

C t {u) = C T (0) + (t - T)(a - A + fa) + o(t - T) 

Now among the four solutions of Remark 3.7, we can choose t and a n such that a 
A + fiii) < 0. In these conditions : 

C t {u) < C T (0) 

which ends the proof. 
2. Final time fixed. 

Definition 3.8. • The trajectory 7 is said to be C° -cost-minimal on [0, T] if there exists a C°-neighborhood 

of 7 such that for any trajectory q contained in this neighborhood, with q(0) = 7(0) and q{T) — j(T), 
we have : Ct{v) > Ct(0), where v is the control associated to q. 

• The trajectory 7 is said to be L°° -cost-minimal on [0, T] if there exists a neighborhood of in L°°([0, T]) 
such that, for any trajectory q associated to a control v of this neighborhood, with q(0) = 7(0) and 
q(T) = 7 (T), we have : C T {v) > C T (0). 

We have the following : 

Theorem 3.9. The trajectory 7 is C° -cost-minimal if and only ifT < t c . Moreover, 7 is not L°° -cost-minimal 
if T > t c (whereas 7 is CP -time-minimal if and only if T < t cc ), where t cc and t c denote the two types of first 
conjugate times ofj. 

Hence in this case, the times at which 7 ceases to be minimizing are different in the time-optimal problem 
and cost-optimal problem : 7 ceases to be C°-time-optimal before it ceases to be C°-cost-optimal (since t cc < t c , 
see Lemma 2.7). 

Proof. The proof is quite similar to proof of Theorem 3.5 (but simpler) and is only sketched. In the coordinates 
(xi, . . . , x n ) of Lemma 2.1, we write E T = (£[, E%), where d 2 E^{0) = 0, and d£f(0), . . ., dE£_i(0), dC T {0) 
are independant. We decompose : 

L 2 ([0, T]) = Rei 8 ■ ■ ■ 8 Me„_i ffite ffiF 

where dEf(0).ei = 1 and dCr(0).eo = 1. By definition of t c : 

• if T < t c : <5/Ke ffiF is positive definite, 

• if T > t c : Q /K eo ©F is indefinite. 

where Q = d 2 E^(0). Hence if T > t c there exist v, w e Ke © F such that Q(v) > and Q(w) < 0. Up to 
adding ee we can assume that projections of v, w on Re are not trivial. Consider the control : 

u = aiei H h a„_ie„_i + Xv + fiw 



- X + f3fi ^ and (t-T)(a- 

□ 
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and try to solve equations : Ei(u) = 8}T, i = 1 . . . n. 

First from the Implicit Function Theorem we get, solving the n — 1 first coordinates (i.e. the controllable 
part of the system) : 

«i = /i(A,/i) = o(A,/i),i = 1.. .n- 1 

Now the last coordinate is : 

E n (u)=Q(u)+o(\\u\\l 2 ) 

= X 2 Q{v) + n 2 Q(w) + X^Q{v, w) + o(A 2 , fi 2 ) 

Using the Mean Value Theorem we obtain : 

• 3Ai < 0, fii > / E n (aiei H h a„_ie„_i + Xiv + mw) = 

• 3X 2 > 0, fi 2 > / E n (aiei H h a„_ie„_i + X 2 v + [i 2 w) = 

There exists i = 1 or 2 such that dCr(0).(Aii> + /U,w) 7^ 0, and therefore there exists a trajectory steering to 
j(T) in time T with a cost strictly lower than Ct(0). □ 

3.2. Optimality of abnormal trajectories for sub-Riemannian systems of rank 2 

With the notations of Section 2.2 we have the following result, which can be proved in the same way as 
Theorem 2.21. 

Theorem 3.10. Under assumptions of Theorem 2.21, the abnormal reference trajectory 7 is (f -optimal for the 
sub-Riemannian system (14) if and only if it is C° -time-minimal for its associated affine system (16). Moreover 
7 is abnormal for this affine system ; actually 7 is C° -optimal if T < t cc and is not L°° -optimal if T > t cc . 

In particular conjugate times are the same along 7 for the sub-Riemannian system (14) and its associated 
affine system (16). Therefore the whole formalism that was introduced for affine systems (the differential 
operators D\,D 2 ) is still valid in sub-Riemannian geometry. Hence the conjugate time of the sub-Riemannian 
problem can be computed using an algorithm. This result makes a link between works of [9] and [4], [5]. 

Example 3.11. The Martinet case (see Section 2.2.2) is in dimension 3, hence t cc = +00 (see Remark 3.3). 
The abnormal trajectory is optimal on R + . 

Remark 3.12. As proved in [3] the C°-optimality is in sub-Riemannian geometry equivalent to the optimality 
in the L 2 topology on controls. 

Remark 3.13. If T is small enough (depending on the choice of the Riemannian structure, and lower than 
t cc ), then as first noted by [4] 7 is moreover globally optimal among all sub-Riemannian trajectories steering 
to 7 (T). 

Remark 3.14. It should be noted that the loss of optimality holds in L°°. Hence using the definitions introduced 
in Section 2.2.2, the L 2 -sector plays no role in the optimality of the reference abnormal trajectory. The loss of 
optimality holds in the L°°-sector. 

I would like to thank my teacher B. Bonnard for many relevant ideas and advices. 
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